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We study numerically the phase-ordering kinetics following a temperature quench of the Ising 
model with single spin flip dynamics on a class of graphs, including geometrical fractals and random 
fractals, such as the percolation cluster. For each structure we discuss the scaling properties and 
compute the dynamical exponents. We show that the exponent a x for the integrated response 
function, at variance with all the other exponents, is independent on temperature and on the presence 
of pinning. This universal charachter suggests a strict relation between a x and the topological 
properties of the networks, in analogy to what observed on regular lattices. 

I. INTRODUCTION 

After quenching a system from an high temperature disordered state to an ordered phase with broken ergodicity, a 
phase-ordering kinetics is observed, characterized by coarsening of quasi-equilibrated domains with a typical size L(t). 
Although a first-principle theory of phase-ordering is presently lacking, for systems defined on homogeneous lattices 
a substantial comprehension of the dynamics has been achieved by means of exact solutions of soluble cases, approx- 
imate theories and numerical simulations While the inner of domains is basically in equilibrium at the quench 
temperature, the off-equilibrium behavior is provided by the slow evolution of their boundaries. As a consequence, 
at each time s, for space separations r <C L(s) or for time separations t — s <C s infra-domains quasi-equilibrium 
properties are probed, while for r 3> L(s) or t — s 3> s one explores inter-domains properties, where the aging behavior 
is observed. Accordingly, the correlation of the order parameter (ai(t)aj(s)) between sites i,j at times s,t can be 
expressed as the sum of two terms 

C ij (t,8)=Cg{t-8)+Cg(t,s). (1) 

The first term describes the equilibrium contribution provided by the interior of domains while the second contains the 
non-equilibrium information. Analogously, also the integrated response function, or zero field cooled magnetization, 
measured on site i at time t after a perturbation has switched on in j from time s onwards, takes an analogous 
addictive form 

Xi j (t,s)=X%(t-s)+X- J 9 (t,s). (2) 

On regular lattices, due to space homogeneity and isotropy, correlation and response function depend only on the 
distance r between i and j. One has, therefore, Cij(t, s) = C(r,t, s), and similarly for Xy (^ s )- 

At the heart of the non-equilibrium behavior is the dynamical scaling symmetry, a self-similarity where time acts 
solely as a length rescaling. When scaling holds, the states sequentially entered by the system are statistically 
equivalent provided lengths are measured in units of the characteristic size Lit) of ordered domains. All the time 
dependence must enter through L(t), and the aging parts in Eqs. (|1I2[) take a scaling form in terms of rescaled 
variables || i = r/L(s) and y = L(t)/L(s) 

C a "(r,t,s)=C(x,y), (3) 

X as (r,t,') = s—*x{x,y). (4) 

The characteristic length grows according to a power law L(t) ~ t x l z . Non-equilibrium exponents such as ci x , z , are 
expected to be universal. Namely, they depend only on a restrict set of parameters, such as space dimensionality and 
number of components J\f of the order parameter, on the conservation laws of the dynamics and are independent on 
temperature. 

In defiance of this basic comprehension in homogeneous systems, our understanding of phase-ordering on inhomo- 
geneous structures is particularly poor, although examples, ranging from disordered materials, to percolation clusters, 
glasses, polymers, and biomolecules, may be abundantly found in physics, economics, chemistry and biology [3[. 
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In this paper we study the phase-ordering kinetics of the Ising model with non-conserved order parameter on some 
physical graphs [4], namely networks with the appropriate topological features to represent real physical structures. 
These structures are constrained to be embeddable in a finite dimensional space and to have bounded coordination 
number. Among these, we consider random fractals, i.e. percolation clusters, and geometrical fractals, such as the 
Sierpinski gasket and carpet and others (see all the cases considered in Fig. Q}. We discuss the results of numerical 
simulations and compare them with the predictions provided by a large- AT model, where exact calculations can be 
carried out [H, @. In the large- M model the general framework of scaling behavior is maintained on generic graphs, 
and the exponents depend on the topology of the network only through the fractal dimension df and the spectral 
dimension d s , a large scale parameter encoding the relevant topology. In particular, one finds z = 2df/d s and 
a x = (d s — 2)/2. This general framework provided by the soluble model seems to be quite unrepresentative of the real 
situation in scalar systems, as observed in 0] and in our simulations. A basic feature missed is the fundamental role 
played by activated processes on inhomogeneous structures. Actually, on homogeneous structures the temperature is 
an irrelevant parameter, in the sense of the renormalization group Long-time, large-scale properties do not depend 
on the quench temperature in the whole low temperature phase. Exponents related to the aging behavior are then 
independent of the thermal noise. On the other hand, structure inhomogeneities exert a pinning force on interfaces or 
other topological defects, whose further evolution can only be mediated by activated processes. Activated processes 
may be present in non disordered models on regular lattices as well, where they sometimes play a fundamental role, but 
without changing large-scale, long-time properties. On fractal sets, on the other hand, pinning forces are exerted on 
all lcngthscales and contribute to the aging behavior. These forces causes a stop- and- go behavior hindering the power 
laws and preventing a straightforward definition of the exponents. At relatively high temperatures, where pinning 
barriers are more easily surmounted, slip-stick effects are less severe. In these cases the general scaling scheme (|3|4[) 
can be investigated and the non equilibrium exponents can be defined. It turns out, however, that pinning forces are 
still subtly at work, making exponents temperature dependent. 

A notable exception in this vague scenario is represented by the response function exponent a x . Interestingly, at 
variance with all the other exponents, its value only depends on the structure considered, is independent on temper- 
ature and on the presence of stick-slip intermittency. This universal character calls for a strict and direct relation 
between a x and precise topological properties of the network, bypassing any microscopic dynamical mechanism. We 
develop an argument showing the relation between this exponent and fundamental equilibrium properties. In partic- 
ular, a x must take a positive value whenever the statistical model on the network considered has a phase transition 
at a finite critical temperature T c , whereas a x = on structures with T c — 0. The same picture is provided by the 
solution of the large- N model on graphs. This is what we find, with good accuracy, in the simulation of discrete 
symmetry models on all the inhomogeneous structures considered. We measure a x = or a x > whenever struc- 
tures with T c = or T c > are considered. Given the fundamental role of topological and connectivity properties 
in determining equilibrium and critical behavior 8] our results provide an evidence for some relationship between 
non-equilibrium kinetics and large scale topology on general networks and suggest that the same topological features 
of graphs determine critical behavior and non-equilibrium exponent a x during phase ordering [f|. 

This paper is organized as follows: In Sec. [IT] we introduce the Ising model that will be considered in the simulations. 
We also introduce the basic observables, and discuss the numerical techniques. In Sec lIIII we discuss the outcome of 
the numerical simulations on different structures. In particular, subsection IIII Bl is devoted to structures with a finite 
T c , namely the percolation cluster above p c , the Sierpinski carpet and toblerone lattices, while subsection IIII Bl deals 
with structures with T c = 0: The percolation cluster above p c , the Sierpinski gasket and the T-fractal. In Sec. IIVI 
we discuss the value of the exponent a x in all the cases considered, and the form of the fluctuation-dissipation plots. 
Sec. [V] contains a final discussion and the conclusions. 



<ij> 

where oi = ±1 is a unitary spin and < ij > are nearest neighbours on a graph. 

The dynamics is introduced by randomly choosing a single spin and updating it with Metropolis transition rate 



II. MODEL AND OBSERVABLES 



The Ising model is defined by the Hamiltonian 




(5) 



»(M -> M) = 5 min [1, exp(-A£/T)] . 



(6) 



Here [a] and [a 1 ] are the spin configurations before and after the move, and 

AE = H[a'} -H[a]. 



(7) 
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Sierpinski Gasket T-Fractal Percolation Cluster 




Sierpinski Carpet Toblerone Lattice Diluted Square Lattice 




FIG. 1: The structures considered in the paper. 

We consider a system of N spins initially prepared in an high temperature uncorrelated state and then quenched, 
at time t — 0, to a low final temperature TV. 

As already discussed, the dynamics of the spins in the bulk of domains provide the equilibrium contribution while 
what is left over accounts for the aging behavior. Since equilibrium dynamics is well understood, the stationary parts 
in Eqs. (I1I2|) are generally well known. In particular, at TV = equilibrium dynamics is frozen and there are no 
stationary contributions. On the other hand, much interest is focused on the aging terms. These can be isolated by 
subtracting the stationary parts computed in equilibrium, from the whole quantity measured in the phase-ordering 
stage. However it is computationally much more efficient to resort to a different method. This amounts to study 
a modified system where TV in the transition rate © is set equal to zero if the spin to be updated belongs to the 
bulk, namely if it is aligned with all its neighbors. Since the bulk degrees of freedom, which alone contribute to the 
stationary parts, feel TV = 0, and equilibrium dynamics is frozen at TV = 0, by computing observables with this 
modified dynamics one isolates the aging term leaving other properties of the dynamics unchanged Q . This modified 
dynamics with no bulk flips (NBF) will be used in the following. The NBF has also the useful property of shifting T c 
towards higher values [10(. As we will discuss below, this fact represents a formidable advantage on inhomogeneous 
graphs. 

Let us discuss the case of graphs where a phase transition occurs at a finite critical temperature T c . We denote these 
as systems of class I. For these structures we can consider quenches to finite temperatures TV < T c . The characteristic 
size L(t) then grows until it becomes comparable with the system size and the new equilibrium state at TV is globally 
attained. We always consider sufficiently large systems to prevent equilibration on the simulated timescales. The 
dynamics in this regime is equivalent to that of an infinite system, for which the final equilibrium state is never 
reached and L(t) keeps growing indefinitely. By analogy with regular lattices, in the late stage one might expect the 
power law 

L(t)~t«, (8) 

although some caveats will be discussed in the following. On regular lattices the value z = 2 is found quite generally. 
In the solution of the large- M model on graphs Q or in the framework of approximate theories Q one finds z = 2df/d s , 
but such a relation is not expected for scalar systems. 

Now we turn to systems of class II, namely those systems for which, with standard dynamics, T c = 0. For every 
finite Tf the system eventually equilibrate to a state with a finite coherence length £(TV), which diverges in the TV — > 
limit. If the temperature is sufficiently low, an interrupted phase-ordering is observed until L(t) becomes comparable 
with £(TV) and equilibration occurs. The phase-ordering phenomenon can be widened at will by decreasing TV, and 
the scaling behavior with Eq. © can be studied. Quite generally, however, on inhomogeneous structures one cannot 
set TV — directly, because this would freeze the dynamics due to pinning effects. For the same reason, also very low 
temperatures are not numerically accessible, since it takes an exponentially long time to surmount pinning barriers. 
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From the numerical point of view, then, one has to find a reasonable compromise between two contrasting issues. 
Namely, Tf must be sufficiently low in order to have a wide scaling regime preceding equilibration, but also sufficiently 
large to allow pinning barriers being overtoped. This may be in some cases impossible. Fortunately, the use of the NBF 
dynamics simplifies considerably the problem. As we have already pointed out, with NBF the critical temperature is 
raised. In the case of system where T c = with standard dynamics, the effect of NBF is to move T c from zero to some 
finite value. Then one has a whole low-temperature region where phase-ordering occurs asymptotically without being 
interrupted by equilibration, much alike in systems of class I. This allows to study the scaling behavior also in this 
case. For all the reasons discussed above, in this paper we will always present results obtained with NBF dynamics. 

Generally speaking, when scaling holds the characteristic length L(t) can be estimated from the knowledge of the 
two-points equal time correlation function, obtained letting t — s in Eq. ([1]) 

G?S(i) = <*(*) •**(*)> (9) 

where (...) means an ensemble average, namely taken over different initial conditions and thermal histories. In 
homogeneous systems according to Eq. © one has 

G%(t) = G a s(r,t)=g(x), (10) 

where x = r/L(t). For scalar order parameter with sharp interfaces, a short distance behavior (x <C 1) of the type 
1 — g(x) ~ x is found [l], [Tlj], namely a Porod's tail g(u) ~ u~( d+1 ) in momentum space for large u — kL(t). With 
Eq. (fTU|) the characteristic length can then be evaluated, for instance, as the half-height with of G(r,i), namely from 
the condition 

G a9 (r = L{t),t) = ^G a9 (0,t) (11) 

On usual lattices the size of domains can be easily related to the density of interfaces pit). In fact, since domains 
are compact, the ratio between their surface and their volume is proportional to L(t) _1 and one has pit) ~ Lit) -1 . 
Since asymptotically pit) has a power law behavior 

P(*)~*-*, (12) 

this implies 

9 = l/z. (13) 

In homogeneous structures, then, p(t) provides an indirect, alternative method for the determination of Lit), and 
hence of z. 

On generic graphs the notion of a distance is not as straightforward as on regular lattices and, considering the 
correlation function © one should, in principle, retain the full dependence of Gi t j(t) on i and j. This would be a 
formidable task without probably providing much insight into the physics and, in particular, into the scaling behavior. 
Therefore, in the following, for some of the structures considered in Sec. IIII1 we will attempt a reasonable definition 
of a distance, we will check the validity of Eq. (jTUJ) with respect to this definition and we will compute Lit) through 
Eq. (flT]) . Let us also notice that the relation (TT3"1) does not hold on generic networks, as we will show explicitly in our 
simulations. Actually, as will be discussed in Sec. IIII A[ the number of spins on the surface of a domain of size Lit) 
may depend dramatically on (among other parameters) the quench temperature. Therefore, in general, there is not 
a unique relation between pit) and L(i) and they provide independent informations. 

The two time quantities that will be considered in this paper are the (spatially averaged) autocorrelation function 

1 N 

and the integrated (auto)response function 

Xit,s) = j dt'R(t,t'). (15) 

J s 

The quantity 



iV 



i ^ aw)) 



dh.it') 



(16) 

h-0 
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is the (spatially averaged) linear response function associated to the perturbation caused by an impulsive magnetic 
field hi switched on at time t' < t. 

On regular lattices, scalings (|3l4p imply 



C a9 (t,s) = h(y) 



and 



X aB (t,s) = s- a *f(y), 
where y = t/s and h(y) = C(Q,y), f(y) = x(0, y), and the large-y behaviors 



(17) 
(18) 
(19) 



f(y) ~ y~ ax 



(20) 



In order to compute \(t, s) we enforce the out of equilibrium generalization of the fluctuation dissipation theorem 
derived in [hS ] , which relates the response function to particular correlation functions of the unperturbed system 



Tx(t,s) 



C(t,t) -C{t,s) 



((Ti(t)Bi(t'))dt' 



where 



Bi[<j] 



(a-i - oQwQcr] -> [a']). 



(21) 



(22) 



In this equation [a] and [a 1 ] are two configurations differing only by the spin on site i, taking the values <7j and a\ 
respectively. This relation allows to compute the integrated response function by measuring correlation functions on 
the unperturbed system, avoiding the complications of the traditional methods where a perturbation is applied, and 
improving significantly the quality of the results [hij . 



III. NUMERICAL SIMULATIONS 



In the following we will present the numerical results. We set J = 1. Statistical errors are comparable to the 
thickness of the symbols. We recall that NBF dynamics is always used, and therefore we measure directly the aging 
part of every observable considered. 



A. Graphs with T c = 



Let start with systems of class II, first. The structures considered in the simulations will be the Sierpinski gasket 
(SG), the T-fractal (TF) and the percolation cluster (PC), namely a diluted square lattice, at the the percolation 
threshold (see Pig. [TJ . The number N of spins in this structures is 265722 and 531442 for the SG, and TS, respectively. 
The PC has been obtained by means of site percolation on a square lattice of size 1200 at the critical dilution p c = 0.407 
The results presented are thermal averages over 100 realizations. 

For the SG we have computed the equal time correlation function ^ by restricting i and j along the borders of 
the structure. This allows a natural definition of the distance r between i and j, since along this lines all sites are 
occupied by spins. For the TF we have computed G ag (r,t) — G"?(t) for points % and j along the the baseline of the 
structure. 

In Fig. [21 G a9 (r,t) is plotted against x = r/L(t) for the SG and the TF. One finds a very good data collapse, as 
expected on the basis of Eq. (fTOjl . This indicates that dynamical scaling is obeyed also on these structures. Notice 
also the presence of the Porod's tail, for small x, implying that interfaces are sharp objects at the relatively low 
temperatures considered in Fig. [21 However, when the temperature is raised, the interfaces broaden on the fractal 
substrate, as shown in Fig. [3] Here one expects a deviation from the Porod linear behavior. This can be clearly seen 
in Fig. [4j Here we plot the exponent c which regulates the small x decay of g(x), namely 1 — g(x) ~ x c , as a function 
of Tf. Clearly, given the choice of i,j described above, our considerations strictly apply only along the particular 
boundary of the structure where G ag {r,t) is computed. However, we expect this results to be representative of the 
whole system, due to the scale invariance of the structure. 
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FIG. 2: G aB (r,t) is plotted against x — r/L(t) for the (TF) and the (SG), the time t ranges from 10 steps for the shortest L(t) 
to 400000 for the longest lengths. 




FIG. 3: (SG) Two spin configurations after t = 1000 step for the temperatures Tf — 0.95 and Tf = 3.45. 



The characteristic size L(i), obtained as in Eq. (fTTj) , is shown in Fig. [5] This length grows roughly as a power law, 
but with a superimposed oscillation which is more pronounced at low temperatures. The origin of such oscillation, 
which is very reminiscent of what one observes in the related problem of Brownian motion on this structure [l5T | , are 
again the pinning forces. In fact, as shown schematically in Fig. [6] for the SG, if two triangles happen to be ordered 
differently, in order to start reversing one of the two to achieve a global ordering, one has to flip one of the interfacial 
spins, for example the one marked with an arrow. This move requires an activation energy AE = 4J and can then 
be accomplished on a time r ~ exp(4 J/Tf). The dynamics then proceeds without activated processes until all the 
triangle is reversed but, then, a new activated step is required and so on. During the time r the dynamics on the 
triangle in consideration is frozen. Since all the triangles have the same qualitative behavior, with some fluctuations, 
the overall behavior of L(t) shows a periodic oscillation due to the recurring slow down caused by activated processes. 
An analogous phenomenon occurs in the TF. Clearly, since r decreases at larger temperatures, this phenomenon is 
less evident, although still present even at relatively high temperatures, as shown on Fig. [5l 

The curves of Fig. [5] can be quite convincingly fitted by a power law with periodic oscillations superimposed. The 
exponent z of the power law, however, should be thought as an effective exponent resulting from the balance between 
non activated processes obeying a growth law ([5]) with a certain value of z, and the activated processes whose net 
effect is an overall slowing down. Since the relative importance of these two processes is regulated by Tf, z turns out 
to be temperature dependent, as shown in Fig. On the basis of the previous discussion, one would expect to have 
a faster growth, namely a larger 1/z, for higher Tf. This is indeed observed in Fig.[7Jin a broad range. For very large 
temperatures 1/z seems to saturate. This is perhaps due to the neighborhood of the critical temperature (that with 
NBF turns out to be T c ~ 3.6 and T c ~ 2 for the SG and the TF, respectively), slowing down the dynamics much in 
the same way as it happens on regular lattices. 

A pattern similar to that of L(t) is observed for pit). This quantity, defined as the density of spin which are not 
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FIG. 4: (SG) The exponent c characterizing the behaviour of g(x) for x ~ 0. 
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FIG. 5: The characteristic length L(t) is plotted against t, for two different temperatures Tf = 1.2 and T/ = 3.0 for the (SC) 
and Tf = 1.25 and T f = 3.0 for the (TF). 




FIG. 6: Schematic representation of a part of the structure, showing the necessity of activated processes. 
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FIG. 7: The dependence on temperature of the exponents 1/z and 8 for the (SG) and the (TF). 
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FIG. 8: The density of interfaces p(t) is plotted against t, at the temperatures T; = 1.2 and Tf = 3.0 for the (SG), at Tf = 1.25 
and Tf = 3.0 for the (TF) and (PC). 



fully aligned with their neighbors, is shown in Fig.[8]for the SG, TF and PC. Also in this case one observes oscillations 
on top of the power law behavior (|12p . The value of 8 also depends on temperature, as shown in Fig. [7] Differently 
from 1/z, however, the behavior is strongly non monotonous, with a broad maximum at intermediate temperatures. 
As already observed, the exponent 9 and 1/z are not trivially related, as on homogeneous structures. In the limit 
of low temperatures one finds an exponent 9 consistent with 9 = —df/z, which can be understood on the basis of 
the following argument. Let us consider again, for simplicity, the SG of Fig. [6l When the temperature is very low 
interfaces are very likely located in the pinning positions, as can be seen in Fig. [3] A domain, therefore, is a triangle 



9 




(SG) 



* s=200 
o s=400 
+ s=1000 
□ s=2000 

* s=4000 





s=200 




s=400 


< 


s=1000 




s=2000 


o 


s=200 




s=400 




s=1000 





s=2000 



(TF) 



T=3 * 





tw100 




tw200 




tw400 


< 


tw1000 





tw100 


o 


tw200 




tw400 




tw1 000 



% 

T=3.00 



(PC) 



FIG. 9: (SG) C a9 (t,s) is plotted against y = t/s, at T) = 3.0 for the (SG), at T f = 3.0 and Tf = 1.0 for the (TF) and at 
T f = 3.00 and Tf = 1.25 for the (PC). 



of size L(t), with a volume V oc L(t) df and a number of surface spins n s 
therefore pit) ~ L(t)~ d f . One then has 



3. The density of interfacial spins is, 



(23) 



instead of Eq. (fT3")) . holding on regular lattices. Let us remark again, however, that, although this relation is consistent 
with our data in the limit of small TV, it is not of general validity. Actually, as already evident from Fig. [3] and from 
the non linear behavior of g(x) at small x, for larger temperatures, interfaces are no longer located on the pinning 
centers. 

Let us consider now the behavior of two time quantities. The autocorrelation function for the SG, TF and PC are 
plotted in Fig. [9] against y = t/s. According to Eq. (fTF)) . one should observe data collapse of the curves with different 
s. This is indeed observed. One also finds the large-?/ power law behavior (119|) with an exponent A strongly dependent 
on temperature. 

Let us turn to consider the response function. According to Eq. (|18|) the exponent a x can be obtained as the slope 
of a double logarithmic plot of x a£, (^ s ) against s, with y held fixed. Such determination should be independent on 
y, within errors. We show this plot in Fig. [10] for the TF. We obtain o Y = .01 ± .02 which is consistent with a x = 0, 
similarly to what found in homogeneous systems with T c = [l|| [TJ [l|[ . Analogous values are found for all the 
structures with T c = 0. Then, by plotting x a9 (£, s) against y one should find data collapse, how it is indeed shown in 
Fig. HUfor SG, TF and PC at different temperatures. 
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FIG. 10: (TF) x a9 (t, s) is plotted against s for fixed values of y = t/s, for T f = 3.00. 
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FIG. 11: (SG) X a9 (t,s) is plotted against y = t/s for different values of s, for T f = 3.00 (SG), T f = 3.00 and 7> = 1.00 (TF) 
and for Tf = 3.00 and Tf = 1.25 (PC). 



B. Graphs with T c > 



In this Section we will consider systems of class I. In particular we will study the phase-ordering kinetics on the 
diluted square lattice (DS) above the percolation threshold, the toblerone lattice (TL) obtained by replicating the 
Sierpinski gasket along a third spatial direction z, and the Sierpinski carpet (SC). These structures are represented 
in Fig. [IJ 

For the SC, we have computed the equal time correlation function along the borders of the structure, analogously 
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FIG. 12: For the (TL) we plotted G% 3 y (r,t) and G a z 9 (r,t) against x = r/L(i) for the (SC) we plotted G a9 (r,t). 



to what done for the SG. The TL is an anisotropic structure, since it is translational invariant along the z direction 
alone. Therefore we have computed the equal time correlation function on the planes with a constant z, and in the 
z direction separately. These are denoted as G"^(r, t) and G" 9 (r, i), respectively. The former is computed along the 
borders of the structure, as for the SG. In Fig. [12] G ag (r, t) is plotted against x = r/L{t) for the TL and the SC. 
One finds a very good data collapse for the TL, as expected on the basis of Eq. (fTU|) . The scaling is less accurate 
for the SG, where the pinning effects seem to be stronger. Notice also in these cases, the presence of the Porod's 
tail, for small x, implying that interfaces are sharp objects at low temperature. As in the case of the SG, when the 
temperature is raised the interfaces broaden and the Porod behavior does not longer hold. 

The characteristic size L(t), obtained as in Eq. (JTTJ) , is shown in Fig. [Jj3 For the TL there are two different lengths 
L xy (t) and L z (t) obtained from G^(r,t) and G^ 9 (r,t) respectively. Note also in this case the presence of oscillations 
on top of an average power-law growth, as discussed in Sec. lIII Al For the SC the exponent z is z ~ 0.33 and z ~ 0.45 
at Tf = 2.00 and Tf = 3.00. The behavior of p(t) is similar to that found for the structures with T c = 0, discussed in 
Sec. UFA] This quantity, is shown in Fig. [H for the DS, TL and SC. For the SC, 9 ~ 0.38 and 9 ~ 0.43 at T f = 2.0 
and Tf = 3.0. Therefore also on these structures both 9 and z depends on the temperature. 

Let us consider now the behavior of two time quantities. The autocorrelation function for the DS, TL and SC are 
plotted in Fig. [15] against y = t/s. 

For the TL one has a good data collapse for the larger values of s, according to Eq. (fT7]) . The situation is different 
for the DS and SC. In these cases, there is no data collapse in the range of time accessed in our simulations. Probably 
this is due to preasymptotic effects and larger values of s then those presented in the figure should be needed to obtain 
the scaling of Eq. (fT7]) . 

Let us turn to consider the response function. From x aff (t, s) we extract a x as the slope of a double logarithmic 
plot of x a9 (t, s) against s, with y held fixed, as described regarding Fig. [TU1 

We obtain a x = 0.21, a x = 0.25 and a x = 0.18 for the DS, the TL and the SC, respectively. Plotting s a * x a9 (i, s) 
against y we find a good data collapse, independently from the temperature, as shown in Fig. [TH indicating that 
Eq. (fT8]) is obeyed. Also the large- y behavior ([20]) is obeyed. 
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FIG. 13: For the TL The characteristic lengths L xy (t) and L z (t) are plotted at the temperatures Tf — 2.0 and Tf = 3.0. For 
the (SC) L(t) is plotted against t, at T/ = 2.0 and T f = 3.0. 



(DS) T=1.25p=0.05 



« T=2 
(TL) a T=3 



* T=2.0 
(SC) o T=3.0 



FIG. 14: (DS) The density of interfaces p(t) is plotted against t, at the temperature Tf — 1.25 for the (DS), at Tf = 2.00 and 
T f = 3.00 for the (TL) and for the (SC). 



IV. THE RESPONSE FUNCTION EXPONENT AND THE FLUCTUATION DISSIPATION PLOT 

In the previous Section, we have checked the validity of the scaling form (fl~8|) in the fractal structures considered, 
and we have determined the value of the exponent a x . As already emphasized, while all the other quantities turn out 
to be very sensitive to the final temperature of the quench, this exponent assumes a well defined value in all the low 
temperature region. On regular lattices, analytical calculations in the vector 0(Af) model in large-A/" limit [19( find 



13 



(DS) 



T=1.25 

O+n * 



s=100 
s=200 
s=400 
s=1000 




(TL) 





s=50 


o 


s=100 




s=200 




s=400 





s=800 




s=1600 




(SO 





s=40 


o 


s=100 




S=200 




s=400 





s=1 000 




s=2000 



FIG. 15: (DS) C ag (t, s) is plotted against y = t/s, at Tf = 1.25 for the (DS) and at Tf = 3.00 for the (TL) and (SC). 



the following dependence on dimensionality 



] d-d L 
du—di 



for 



d < d L 



with log corrections for 
for d > djj, 



d = dr 



(24) 



where d^ is the lower critical dimension of static critical phenomena and djj = 4. Numerical simulations [13, l2Q|, l2l(| of 
discrete and 0(M) vectorial systems, with conserved and non conserved order parameter, are consistent with Eq. (|24p. 
and djj = 3 or djj — 4 for scalar or vectorial models respectively. This shows that the non-equilibrium exponent, 
a x , is related in a non-trivial way to the topology of the underlying lattice, through d. Interestingly, this relation 
implies a x > when the system is above c?l, i.e. when a phase transition at finite temperature T c is present, a x = 
at dL and a x < below d^. This result for regular lattices suggests that the non-equilibrium dynamics of statistical 
models, and in particular the response function exponent, could be related to important topological properties also 
in the case of generic graphs. This hypothesis can be checked in some detail in systems with a continuous symmetry, 
because in this case it is well known that a large scale parameter, the "spectral dimension" d s , encodes the relevant 
topological features, uniquely determines the existence of phase transitions [22j and controls the critical behavior [23| . 
In other words on generic networks d s plays the same role played by the Euclidean dimension d on regular lattices. In 
particular, vectorial models on graphs exhibit a phase transition at finite temperature only if d s > dL = 2. Solving an 
0(J\f) model in the large- Af limit [6| one can show that this property holds true also for the non-equilibrium exponent 
a x . In fact, one finds the same expression (|24[) as for regular lattices, with d s replacing d. Then, again one has a x > 
when the system is above di, i.e. when a phase transition at finite temperature T c is present, and a x < when T c = 
The conclusion is that, for models with a continuous symmetry there is a well defined relation between the non 
equilibrium exponent a x and the topological features of the network which regulate the critical properties. In this 
case, this features are totally described by a single index, namely d s . The next question is if a similar picture hold 
for scalar systems, where the counterpart of d s is not known. Namely, if the non equilibrium exponent a x is related 
to the critical properties, and in turn to topology, such that a x > or a x < holds for graphs with T c > or T c = 
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FIG. 16: X aa (t, s) is plotted against y = t/s for different values of s, at T> = 1.25 for the (DS), at T; = 2.00 and T f = 3.00 for 
the (TL) and the (SC). 



respectively. An argument based on the topology of class II networks suggests that a x = is expected for these 



systems. The values of a x reported in the previous Section IIIII show quite convincingly that the above conjecture is 
verified in all the cases considered in our simulations. One has a x ~ in all the structures of class II and a x > in 
those of class I. This result is particularly interesting for discrete symmetry models. In fact, while for 0(JV) models 
the presence of phase transition can be inferred from d s , there is not a such a general criterion for discrete models. 
Our data suggest that a x may be used to infer the presence of phase transition on a generic network. We recall again 
that this result directly links a x to some relevant topological features of graphs. For this reason, although a x is a 
non-equilibrium dynamical exponent, it is totally independent on all the non-universal details of the dynamics we 
have described in Sec. IIIII We have already observed, in fact, that, differently from all the other dynamical exponent, 
its value is robust and temperature independent. 

Finally, we mention that the value a x = found in the structures of class II is associated to an anomalous 
fluctuation-dissipation plot [2(|. Re-parameterizing in x{t, s) the time t in terms of C(t, s), according to the theorem 
by Franz, Mezard, Parisi and Peliti [2J|, under broad hypoteses the following relation 



Tf lim 

S — >OG 



d 2 x(C, s) 



dC 2 



(25) 



C=q 



exists between the non equilibrium two time functions x, C and the equilibrium probability distribution P eq (q) of the 



overlaps Q([cr], [a 1 ]) = (1/N) J2i=i a i a 'i between two configurations [a], [a'] 

P eM = i E exp{-W[*] + H[a'])}S(Q([al [a']) - q). 



(26) 



Using the 8-\ike form of P e q(q) of the low temperature state of the ferromagnetic models considered in this paper, one 
obtains 20] the well known broken line fluctuation-dissipation plot of x versus C. However, as discussed in [6J, when 
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a x = the theorem (|25|) cannot be applied and one gets a non-trivial fluctuation-dissipation plot which is not related 

to P eq (q)- 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have studied the phase-ordering kinetics of the Ising model with spin flip dynamics on a class of 
physical graphs. The evolution is in many respects qualitatively similar to what observed on regular lattices. After the 
quench one has the formation and growth of domains of the two competing ordered phases, and a scaling symmetry is 
obeyed quite generally (except for the SC in the timescale of our simulations). However, differently from coarsening 
on regular lattices, the fractal nature of the networks pins the interfaces on locally stable positions. Escape from the 
pinned positions is achieved by means of activated processes, but subsequently the interfaces are trapped again. When 
interfaces are pinned, the growth of the domain size L(t) is inhibited. Because of this recurrent phenomenon the usual 
power growth-law of L(t) is modulated by an oscillation, which is more evident when the pinning is stronger, namely 
at low temperatures. The necessity of activated moves on all time/space scales makes a great difference with respect to 
regular lattices. On regular lattices the temperature is an irrelevant parameter [l[, in the sense of the renormalization 
group. Namely, in all the low temperature region the dynamics is controlled by the zero temperature fixed point. One 
has, therefore, the same dynamical exponents in all quenches to Tf < T c , regardless of the temperature. On fractal 
networks, instead, the situation is different, because the dynamics of activated processes are faster the higher is Tf. 
As a consequence, exponents are always found temperature dependent. 

A notable exception is the response function exponent a x , which turns out to be stable and temperature independent. 
This already suggests that a x may be related to a more fundamental property of the system which remains stable 
under temperature variations. Interestingly, we find a x — on all the considered graphs which do not support a 
ferro-paramagnetic transition at a finite critical temperature, while a x > in all the other cases. The same situation 
was found [H, [2(| HH on regular lattices, where a x is related to the euclidean dimension. On these lattices, the 
euclidean dimension is the topological parameter that determines the exhistence of phase transitions and regulates 
the critical properties. This calls for the hypoteses that, also on generic graphs, a x could be related to the relevant 
topological features which govern the critical behavior, although for systems with a discrete symmetry, such as the 
Ising model considered here, a unique topological indicator, analogous to d is not known. Interestingly, this fact 
suggests that a x can be used to infer topological properties of graphs. 
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